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ABSTRACT 

We study the physical properties, formation histories, and environment of galaxies without 
a significant "classical" spheroidal component, as predicted by semi-analytical models of 
galaxy formation and evolution. This work is complementary to the analysis presented in 
iDe Lucia et al.l (1201 Ih . where we focus on the relative contribution of various physical mech- 
anisms responsible for bulge assembly in a ACDM cosmology. We find that the fraction of 
bulgeless galaxies is a strong decreasing function of stellar mass: they represent a negligible 
fraction of the galaxy population with M^, > lO^^A/©, but dominate at M^, < 10^"^/©. We 
find a clear dichotomy in this galaxy population, between central galaxies of low-mass dark 
matter haloes, and satellite galaxies in massive groups/clusters. We show that bulgeless galax- 
ies are relatively young systems, that assemble most of their mass at low-redshift, but they 
can also host very old stellar populations. Since galaxy-galaxy mergers are assumed to lead 
to the formation of a spheroidal component, in our models these galaxies form preferentially 
in low-mass haloes that host a small number of satellites galaxies. We show that the adopted 
modelling for galaxy mergers represents a key ingredient in determining the actual number of 
bulgeless galaxies. Our results show that these galaxies are not a rare population in theoretical 
models: at z ~ 0, galaxies with no classical bulge (but often including galaxies with the equiv- 
alent of pseudo-bulges) account for up to 14% of the galaxies with 10^^ < A/^/A/© < 10^^. 

Key words: galaxies: formation - galaxies: evolution - galaxies:bulges - galaxies:interactions 
- galaxies:structure 



1 INTRODUCTION 

Sin ce the in t roduc tion of the morphological classification scheme 
by iHubbld jl92d) . two components are traditionally identified 
in galaxies: a centrally concentrated spheroidal-like structure 
("bulge") and a disc-like stellar distribution, often associated with 
spiral arms. The galaxy population can be (and usually is) clas- 
sified according to the relative contribution of these two compo- 
nents to the total light emitted by the system. A finer classification 
takes into account the contribution from other components (e.g. 
bars, spiral arms). In the last few decades, observational evidence 
has been gathered to indicate that this picture is oversimplified: 
bulges (which contribute up to 60 % of the stel lar mass in mas- 
sive galaxies in the local Universe, lGadottiir2009h are now seen as 
a heterogeneous class, including purely spheroidal systems (ellipti- 
cal galaxies), "classical" bulges (dynamically and photometrically 
similar to ellipticals, but with sign ificant kinematical differences, 
see e.g. lDavies &Illingworti]| 19831) , and "pseudo" bulges (charac- 
terised by "disc"-like expone ntial profiles or kinematics, see e.g. 
iKormendv & Kennicutl|[2004l and references herein). 



Theoretical models predict that early star formation takes 
place mainly in discs that form due to the conservation of the an- 
gular momentum acquired through early torques acting during the 
proto-galactic stage. Bulges form as the result of physical processes 
able to remove angular momentum from stars and gas. In partic- 
ular, classical bulges and purely spheroidal systems are believed 
to be associated with the most violent dissipative processes, like 
mergers and close interactions. On the other hand, pseudo bulges 
are usually linked to secular evolution of gravitational instabilities 
in the disc component: the "unstable" structure is expected to find 
a new equilibrium following a rearrangement of part of the gas and 
stars in a central structure with enhanced density. Whenever these 
processes are infrequent and/or inefficient, we expect the galaxy 
morphology to be dominated by a disc component. These galaxies 
are often referred to as "bulgeless", and have been seen as a poten- 
tial challenge for cuiTent theor ies of galaxy formation and evolution 
in a ACDM Universe (see e.g. lD'Onghia & Burkertl2004l). Follow- 
ing th ese suggestions several authors (see e.g. Graham & Worle^j 
l2008h used bulge-disc decomposition algorithms to determine the 
relative contribution of these two components in samples of local 
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galaxies. In particular. lWeinzirl et al.l ( l2009h considered a sample of 
143 bright {Mb < —19.3) low inclination spirals and found that 
a relevant fraction (~ 66%) of these is dominated by a disc com- 
ponent accounting for more than 80% of their total stellar mass. 
Starting from a sam ple of ~ 4000 bright ( Mb < —18) galaxies at 
0.013 < z < 0.18. lcameronetal.l ( l2009h computed a total stellar 
mass de nsity for pure-dis c objec ts of 1.3 ± 0.1MqMpc~^ . More 
recently, iKormendv e7ai] l l201(]l KIO hereafter) consider a sample 
of 19 relatively massive (rotation velocities Kot > 150km/s) and 
close (< 8 Mpc) galaxies, and find that 4 of these are consistent 
with being pure disc galaxies. 7 galaxies in the same sample (in- 
cluding the Milky Way) have pseudo-bulges. In total, KIO estimate 
that 58 — 74% of the galaxies in their sample do not experience 
relevant mergers in the past (they include in the bulgeless cate- 
gory those galaxies with a pseudo-bulge). Considering that, at these 
masses, quiet merger histories are rare, all these authors propose 
their estimates as a crucial test for galaxy formation models. 

The formation of bulgeless galaxies represents a classical 
challenge also for cosmological N-body hydrodynamical simula- 
tions. If conserved during its collapse toward the centre of the 
galaxy, the angular momentum of the gas is sufficient to produce 
large discs (e.g., fall & Efstathiou 1980; Mo et al. 199^. Cool- 
ing in small and dense progenitors at high redshift, however, con- 
dense the gas in their inner regions. Dynamical friction on the or- 
biting satellites then dissipates the gas angular momentum (e.g., 
iD'Onghia et"ai]|2006l) . As a result, discs in N-body simulations are 
too compact (with rotation curves that usual ly peak at a few kpc) 
with resp ect to observationa l measurements dSteinmetz & Navarrol 
1 19991 : see lMaver et alj2008l for a review). One of the proposed (and 
so far most successful) solution to the angular momentum 'catas- 
trophe' requires a fine tuning of the feedback in s mall high-redshift 
proge nitors so as to avoid early gas cooling (e.g.. lGovemato et al.l 
I2OI0I) . It should be noted that loss of angular momentum is caused 
also by secular evolu tion of discs and bar instability (see e.g. 
lOebattista et al.|[2006l) . and that these are more easily triggered in 
discs that are mad e compact by partial loss of angular momentum 
dCurir et al.ll2006l) . 

In a recent paper, lOe Lucia et al.l ilOl iL hereafter Paper I) 
analyse the predictions from semi-analytical models (SAMs) of 
galaxy formation and evolution within a ACDM cosmology, and 
quantify the relative contribution of different processes (major and 
minor mergers, and disc instabilities) to the assembly of bulges. In 
this paper, we will tackle a complementary question, i.e. determine 
under which conditions a model galaxy does not develop a signifi- 
cant spheroidal component. 



2 MODELS 

In this paper, we present p redictions from two independently de- 
veloped SAMs, namely the lWang et al] j2008l , hereafter WDL08) 
implementation of the "Munich" model and the MORGANA m odel, 
as adapted to the WMAP3 cosmology in ILo Faroetal] ( l2009l) . We 
refer to the original papers for a detailed discussion of the mod- 
els used in this study, and to Paper I for a detailed description of 
the recipes adopted to model bulge formation. In this section, we 
provide a brief summary of these ingredients. 

Both models consider similar channels leading to the assem- 
bly of a spheroidal component, namely galaxy-galaxy mergers and 
disc instabilities. There are however, significant differences in the 
treatment of these processes. In both models, major mergers com- 
pletely destroy the disc components of the two galaxies. The rem- 



nant spheroidal galaxy may eventually regrow a disc, if fed by 
an appreciable cooling flow at later times. During minor merg- 
ers, both SAMs assign the stellar component of the secondary 
galaxy to the bulge of the remnant, but make different assumptions 
about the stars formed during the burst associated with the merger: 
WDL08 gives them to the disc of the remnant galaxy, while MOR- 
GANA gives them to its bulge. This different choice implies that 
the contribution of mi nor mergers to bulge assembly is more im- 
portant in MORGANA. iDe Lucia et al.l 1 120 id) compared the merger 
models implemented in MORGANA and WDL08 and found that the 
former provides merger times that are systematically shorter than 
those used in the latter (by about an order of magnitude): this trans- 
lates into a higher frequency of merger events in MORGANA and 
has important consequences for the timing of bulge formation, as 
discussed in Paper I. In particular, due to the shorter merger times 
and the different treatment of minor mergers, we have shown that 
MORGANA predicts a larger stellar mass locked in bulges at each 
redshift, and larger mean bulge-to-total ratios for galaxies of all 
masses. 

The two models used in this study also differ in their treatment 
of disc instabilities: both models adopt the stability criterion de- 
fined in lsfstath iou et al. (1982), but use different definitions for the 
relevant physical quantities (in particular disc velocities and radii, 
see sec. 7 and Fig. 10 of Paper I). In addition, they make different 
assumptions about the re-arrangement of baryons following insta- 
bility events: in WDL08 only the stellar mass fraction necessary to 
restore stability is transferred from the disk to the bulge. In MOR- 
GANA, a significant fraction (i.e. half) of the baryonic mass (both 
gas and stars) of the disk is transferred to the bulge. As shown in 
Paper I, the approach adopted in the MORGANA model translates 
into a more prominent role of the disk instability channel in bulge 
formation. 

Following Paper I, we consider in the following three differ- 
ent implementations for each model: a standard implementation, 
which includes both mergers and disc instability, a pure merger 
implementation where we sw itch off the disk insta bility channel, 
and a model that adopts the iHopkins et"an j2009l HOP09 here- 
after) prescriptions for the re-distribution of gas and stars during 
mergers, and for modelling the fraction of disc material which sur- 
vives merger events. Briefly, the HOP09 approach reduces the ef- 
ficiency of bulge formation and increases the fraction of baryonic 
mass in disc components at each redshift. In addition, it assumes 
that a fraction of the disc survives even during major mergers. In 
Paper I, we showed that this change has important consequences for 
the predicted space density of purely spheroidal (elliptical) galax- 
ies, but it does not affect bulge formation in galaxies less massive 
than~ 10" Mq. 

Our SAMs do not allow a fine classification of the different 
bulge subclasses to be made. In particular, it is not possible to disen- 
tangle between the formation of a pseudo or classical bulge, based 
just on the properties of the final spheroidal component. Nonethe- 
less, our strategy provides a natural framework for the analysis 
presented in this study: assuming that classical bulges are associ- 
ated with mergers, and that pseudo-bulges originate from instabil- 
ities, the standard implementation gives the full statistics for disc- 
dominated galaxies, while the comparison with the pure merger im- 
plementation provides information about the relative contribution 
of classical and pseudo-bulges. 

In this paper we will define as "bulgeless" all model galax- 
ies with a bulge-to-total (B /T) mass ratio lower than 0.1, i.e. all 
model galaxies whose bulges contribute to less than 10% of their 
total stellar mass. In the following, we consider only galaxies with 
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Figure 1. The fraction of bulgeless galaxies (/bi) as a function of their stel- 
lar mass and parent halo mass. In each panel, red, blue and green lines refer 
to the standard model, to the pure merger implementation, and the HOP09 
implementation respectively. Shaded histograms are for central galaxies, 
and use the same colour coding. In the panels showing results from MOR- 
GANA, the thick solid black line shows results obtained adopting longer 
merger times (see text for more details). 

Mi, > IQ^ Mq, which is above the resolution limit of our sim- 
ulations. Resolution effects may play a role in our definition of 
bulgeless galaxies, especially for galaxies with M* < IO^'^A/q, 
since we are not able to account for mergers with and between 
lower mass galaxies. We note, however, that galaxy formation in 
lower masses haloes is an increasingly inefficient process (see e.g. 
ISenson & Devereu^l2010l) . So minor mergers with small galaxies 
hosted in sub-resolution haloes are not expected to contribute much 
to the growth of the central spheroid, and to affect the results dis- 
cussed below. 



3 RESULTS 

In this section, we will discuss in detail the properties, environ- 
ment and formation history of bulgeless galaxies, as predicted by 
our SAMs. In particular, we consider the fraction fbi{P) of model 
galaxies with B/T < 0.1 and with a given property P, and the 
normalised (to the total number of systems) distribution nbi(P) of 
bulgeless galaxies as a function of P. 

Fig. [T] shows /bi(A/) as a function of stellar mass (Af*) and 
parent halo mass (A/dm). In each panel, the red, blue and green 
histograms refer to the standard, pure merger and HOP09 imple- 
mentations respectively. The shaded histograms show the contri- 
bution of central galaxies (same colour coding). Common trends 
between the two SAMs and the three implementations considered 
can be seen. In the standard runs, fhi{M-t) exhibits a marked de- 
crease as a function of Af* (left panels). A similar result is obtained 
in the pure merger runs, although in this case the resulting fractions 
are larger than those obtained in the standard models, because of 
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Standard Implementation 
WDL08 26l" % 


otl" % 


KIO (total bulges) 60^33 % 




Pure Mergers Implementation 
WDL08 57l" % 


utll % 


KIO (classical bulges) SO^gg % 


55^23 % 



Table 1. Predicted and observed fractions /bi(Aft) of bulgeless galax- 
ies around MW-like haloes. Theoretical predictions refer to the median 
/bi{Afi) value, with the confidence levels defined on the fith and 95*'' 
precentiles of the distribution, while confidence levels for the KIO data are 
based on the Wilson score interval approximation. 



the lack of bulges forming via the disc instability channel. Neglect- 
ing bulge formation via disc instabilities has a stronger effect on 
predictions from MORGANA than on those from WDL08. This is 
due to the more efficient mass transfer associated with this phys- 
ical mechanism in the former model. At all masses, /bi(Af*) in 
WDL08 is larger than in MORGANA for both implementations: this 
is due to the assumption of shorter merger times in MORGANA, 
that leads to a more efficient bulge formation in this model with 
respect to WDL08. Therefore, the relatively small /bi(A/*) in the 
standard MORGANA run is due to a combination of shorter merger 
time scales, and stronger mass transfer via disc instabilities. We 
test explicitly the influence of different merger times, by re-running 
MORGANA using the same dynamical friction prescription adopted 
in WDL08 - results are shown as a black solid line in the bottom 
panels of fig. [T] As expected, when longer merger times are as- 
sumed, /bi(Af) increases because more galaxies are able to avoid 
mergers. 

We find that the HOP09 implementations do not alter signif- 
icantly the fraction of bulgeless galaxies with respect to the runs 
considered above: this is consistent with conclusions from our Pa- 
per I that this recipe does not affect bulge formation in galaxies less 
massive than ~ IO^^A/q. We note, however, that the HOP09 pre- 
dictions for the WDL08 model are closer to the standard results, 
while for MORGANA they are closer to the pure merger model re- 
sults. This difference is due to the different treatment of the stars 
associated with bursts triggered by minor mergers: as explained 
above, these are given to the disc of the remnant galaxy in the 
standard implementation of the WDL08 model, while when adopt- 
ing the HOP09 recipes this model assumes (as done in the MOR- 
GANA model) that these stars go to the bulge of the remnant galaxy. 

We then consider the fraction of bulgeless galaxies as a func- 
tion of their parent halo mass (/bi(A/DM), right panels). A clear 
dichotomy in the bulgeless population can be seen: lower-mass 
haloes have an increasing probability of hosting a central bulge- 
less galaxy (larger than 70% for A'/dm < 10^^ Af©), while almost 
all central galaxies of haloes with A/dm XL 10^^ Mq host signif- 
icant bulges in all implementations. On the other hand, bulgeless 
satellites constitute an important fraction of cluster/group galax- 
ies. MORGANA predicts a lower /bi(AfDM), and a clear dip at 
A/dm ~ lO^^A/0 (depending on the chosen implementation). For 
the WDL08 model the trends as a function of the halo mass are 
somewhat weaker. 

In order to compare the results shown in fig. [T]with the KIO 
data, we have computed the galaxy stellar mass for all galax- 
ies in the KIO sample. In particular, we have used table Bl 
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from IZibetti et alj ( |2009|) . the A'-band magnitudes given in KIO, 
and (B — V)(i col ours obtained using the HyperLeda database 
dPaturel et al ] i2003h . Given the tight relation between the stellar 
mass of central galaxies and their parent halo mass, we expect 
this approach to be consistent with the KIO analysis (which is 
based on rotation velocities), and it provides a more straightfor- 
ward comparison with our theoretical predictions. We thus com- 
pute the fractions of bulgeless galaxies /bi(Af*) in two mass bins 
(10^° < M^/Mq < 10", and lO" < M^,/Mq < 10^^), by 
using the B/T ratios listed in KIO (see their table 2). We give our 
estimated fractions in Table[T] together with the Wilson score inter- 
val|3 calibrated at the 95% confidence limit. 

In order to take into account the error due to cosmic variance, 
we consider all central galaxies living in Milky Way-type haloes 
^Qi2±o.25^^^ in the WDL08 model, and define KlO-like sam- 
ples in the mass range lO'^" < Mum/Mq < 10^^ by consider- 
ing all galaxies closer than 8 Mpc. We then compute the median 
/bi(Af*) value and its confidence interval based on the fifth and 
95**^ percentiles of the distribution. This analysis is limited to the 
WDL08 model because MORGANA does not predict accurate po- 
sitions for satellite galaxies. In Table [T] we compare predictions 
for the standard WDL08 run with the fraction of disc galaxies in 
KIO sample (i.e. galaxies with total B/T < 0.1), and predictions 
obtained from the pure-merger model with the fraction of galaxies 
with no classical bulge defined in KIO (i.e. classical B/T < 0.1). 
The theoretical fractions are systematically lower than observa- 
tions, and the discrepancy is severe for the standard implementa- 
tion compared to the distribution of pure discs: the probability of 
observing a sample with the same morphological mix measured by 
KIO is smaller than 1% in the volume of the simulation used by the 
WDL08 model. The discrepancy between theoretical predictions 
and observational estimates is, however, reduced if we consider the 
predictions for the pure merger runs and compare them to the frac- 
tion of galaxies without a classical bulge: in this case, the probabil- 
ity of finding a MW-like neighbourhood similar to the KIO sample 
is larger than 5%. We note that when selecting MW-like haloes, 
we have not applied any isolation criterion, which could further in- 
crease the expected fraction of bulgeless systems. If pseudo-bulges 
can be associated with secular processes, our results clearly show 
the importance of the adopted modelling of disc instabilities in or- 
der to correctly estimate the expected /bi(M). 

We also compare pre dictions from theoreti cal models with the 
observational estimates o f lWeinzirl et aLl(l2009h and lCameron et all 
( l2009h who compi led samples that are larger than KIO. How- 
ever, the iWeinzir l etal] ( |2009|) sample is not mass and volume 
complete due to less well controlled selection criteria, while the 
ICameron et al.l ( l2009h decompositions suffer a number of addi- 
tional problems, due to poorer resolution of higher redshift objects: 
in fact we note that accurate decomposition of galaxies is difficult, 
becoming worse at low resolution (e.g . higher z) and where fe wer 
components are included (e.g. bars, see lLaurikainen et al .120071) . In 
addition, an accurate comparison with these samples would require 
additional sources of uncertainties in our model predictions, such 
as the inclusion of dust attenuation and the modelling of synthetic 
SEDs. 

Nonetheless, we make a qualitative comparison between data 
and models, by defining a mass selected sample of Af* > 10^° A/© 



1 The lWilsoj 1 1921) approximation provides an estimate for asymmetric 
errors based on confidence levels of a binomial distribution, which is more 
appropriate than a normal distribution for small number statistics. 
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Figure 2. The fraction /bi of bulg eless galaxies a s a fun ction of stellar 
mass. Observational constrains from lFisher & Drorvl fcoill) data are shown 
as solid (total B/T < 0.1) and dashed lines (classical B/T < 0.1), while 
Red and blue dots refer to the standard and pure merger implementations of 
the WDL08 model (see text for more details). 

model galaxies, which roughly corresponds to the mass range cov- 
ered by observations. We find the following z = total stel- 
lar mass densities for B/T < 0.1 galaxies: 0.3 and 0.9 x 
l{f MqMpc~^ (for the standard and pure merger implementa- 
tions of the WDL08 model, respectively) and 0.03 and 0.6 x 
lO^A^Q A'/pc^^ (for the standard and pure merger implementations 
of MORGANA model, respectively). T hese numbers are lower than 
those found bv lCameron e7aLl ( l2009l 1.3±0.1 x l(f M^Mpc-^). 

Using the same sample of model galaxies, we compare the 
predicted frac t ion o f B/T < 0.2 high-mass spirals with the 
IWeinzirl et all ( |2009|) results (~ 66%, following that paper, we de- 
fine spirals as galaxies with B/T < 0.75). These galaxies account 
for 18% (67%) of objects in the standard (pure merger) implemen- 
tation of the WDL08 model, and for 49% (81%) of objects in the 
standard (pure merger) implementation of MORGANA. These esti- 
mates confirm our conclusion that bulgeless galaxies are underpre- 
dicted by the standard implementations of semi-analytical models. 
The discrepancy is, however, sensibly reduced when pure merger 
implementations of the same models are cons i dered. 

In a recent publication, Ipisher & Drorvl (|20I lb analysed the 
fraction of classical and pseudo-bulges in the local volume (closer 
than II Mpc). We repeat our analysis on /bi(Af*), using this new 
sa mple and the same appro ach as for the KIO sample (fig.|2]l. From 
the lFisher & Drorvl 1 I2OI ij) data (their table. 1), we draw two sub- 
samples including galaxies with both total B/T < 0.1 (solid line 
and gray shaded region) and classical B/T < 0.1 (dashed line, 
we consider all galaxies showing a prominent pseudo-bulge com- 
ponent as bulgeless), and we compare them with the prediction of 
the standard (red diamonds with errorbars) and pure merger (blue 
diamonds with errorbars) implementations of the WDL08 model, 
respectively. Samples of model galaxies are defined within 1 1 Mpc 
from central galaxies of Milky Way-like haloes and we compute 
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Figure 3. Upper panels: /bi as a function of the number of effective satel- 
lites (see text for more details). Lower panels: n^i as a function of the num- 
ber of effective satellites. In each panel colours, lines and shadings are the 
same as in fig.[T] 



the mean /bi(M*) and its confidence interval based on the S**" and 
OS**" percentiles of the distribution. In both cases, model predic- 
tions underpredict the observational results; however, the discrep- 
ancy between the standard implementation and the /bi(M*) based 
on the total B/T < 0.1 ratio is reduced with respect to the KIO 
sample. The agreement between the pure merger implementation 
and the sample based on the classical B/T < 0.1 is satisfactory. 
We note that theoretical predictions obtained when considering the 
entire simulation box are very clos e to the mean shown i n fig. |2] 
It is also worth mentioning that the lFisher & Drorvl ( 1201 ll) sample 
contains a significant number of galaxies hosting a pseudo-bulge 
and with total B/T > 0.3 (they account for ~ 10% of the whole 
sample of galaxies hosting a pseudo-bulge). Theoretically, we ex- 
pect that such prominent pseudo-bulges may hide a sub-dominant 
(hence relevant) merger-driven classical bulge component and this 
effect would reduce the discrepancy between the dashed line and 
blue diamonds. 

In order to provide limits on the number of mergers suffered 
by bulgeless galaxies, we take advantage of the merger histories 
provided by our models. In particular, for each galaxy in the model, 
we define a number of "effective satellites" (A'sat) as the number of 
satellite galaxies in its parent halo, computed at the redshift when 
it (its main progenitor) was for the last time a central galaxy. We 
only consider satellite galaxies with A/* > IO'^'Mq. For central 
galaxies, the number of effective satellites corresponds to the num- 
ber of actual satellites. We show /bi( A'sat) in fig.[3](upper panels). 
We note that the probability of being bulgeless increases with de- 
creasing TVsat- The Figure shows that MORGANA exhibits a very 
skewed distribution towards low A'sat values, while the con'espond- 
ing distribution for WDL08 is flatter. This is a consequence of the 
different merger times adopted in the two models, as confirmed by 
the solid thick line in the right panels that coiTespond to a run of 



the MORGANA model (the pure merger run) with longer merger 
times. Since MORGANA assumes (in its standard implementation) 
quite short merger times, only galaxies living in very isolated envi- 
ronments survive as bulgeless systems. The WDL08 model adopts 
longer merger times, so that also galaxies living in haloes with more 
substructures may avoid developing a significant bulge. Disc insta- 
bilities do not change the distributions, but as expected, they affect 
the fraction of bulgeless galaxies. We stress that, despite a wider 
range of allowed halo merger histories for bulgeless galaxies in 
WDL08, the distribution of nbi (A'sat) (lower panels) is peaked to- 
wards very small A'sat values. 

In order to better characterize the physical properties of bulge- 
less galaxies as predicted from our models, we consider the quan- 
tity Zhaif, defined as the redshift at which half of the final stellar 
mass is assembled in a single object. Since mergers play a small 
role in the assembly of bulgeless galaxies, this quantity is also a 
good indicator for the star formation history of the galaxy (i.e. 
is a good approximation of the time when half of the stars were 
formed). The distribution of 2haif predicted by our models is shown 
in fig.|4](left panels), and shows that bulgeless model galaxies cor- 
respond to a relatively young population: central galaxies are more 
skewed towards low Zhaif with respect to satellites. These results 
are consistent with recent observations of the central regions of 
M33, whi ch appear to be d ominated by an old stellar population 
6 Gvr |javadiet"Zll201 ll) It is also worth noting that these ob- 
jects contain a significant fraction of old stars in their discs: the 
presence of 9-10 Gyr old stars in a Milky Way-like galaxy is not 
unusual. Removing disc instabilities or using the HOP09 merger 
prescription, does not modify significantly the Zhaif distribution: 
this is due to the paucity of bulge forming events in the bulgeless 
galaxies' history. We also consider the distribution of accretion red- 
shifts for the bulgeless satellites (zsat, fig.|4]right panels): here we 
define Zsat as the last time the galaxy is the central object of an 
independent DM halo. It is worth stressing that this definition does 
not always correspond to the redshift when the satellite is accreted 
onto the main progenitor of its 2 = parent halo. Nevertheless, 
since Zsat the star formation history of satellite galaxies is strongly 
affected by the strangulation in both models. The Figure shows that 
the overall Zsat distribution is broader than the Zhaif one. There- 
fore, both models predicts a population of pure red discs in groups 
and clusters, dominated by quite old stellar populations, which have 
been accreted as satellites and have stopped forming stars relatively 
recently. 



4 DISCUSSION & CONCLUSIONS 

In this paper we study the statistics, distributions and formation 
histories of galaxies with bulge-to- total mass ratios B/T < 0.1, 
as predicted by theoretical models in the ACDM f ramework. This 
is the second paper of a series: in the first paper dOe Lucia et al.l 
llOllI) . we have studied the formation of spheroids, and in a third 
paper we will focus on a detailed comparison between model pre- 
dictions and observational data (Wilman et al, in preparation). 

We conside r two independently developed SAMs: the 
IWang et al] ( I2OO8I) implementation of the "M unich" model, and th e 
most recent implementation of MORGANA jMonaco et al.ll2007h . 
In addition, we consider three different implementations of each 
model: a standard run, a pure merger model, and a modified merger 
model based on the prescriptions proposed bv lHopkins et al. I( l2009l) 
and based on recend hydrodynamical simulations of galaxy merg- 
ers. The two models used in this study include the same channels 
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Figure 4. Left panel: the distribution of formation redshifts of bulgeless 
galaxies (see text for more details). Right panel: tlie distribution of accretion 
redshifts for bulgeless satellites. In each panel colours, lines and shadings 
are the same as in fig.[T] 



for bulge formation (i.e. galaxy-galaxy mergers and disc instabili- 
ties). 

Our results highlight that models predict a non-negligible frac- 
tion of bulgeless Af* < 10^^ galaxies at 2 ~ 0. For all models 
and implementations considered, the fraction of bulgeless galaxies 
decreases rapidly with increasing stellar mass, and becomes neg- 
ligible for M* > 10^^ Mq. Bulgeless galaxies are either central 
galaxies in low mass haloes, or satellites in groups and clusters. 
They assemble their mass at relatively low redshifts, but they can 
host quite old stellar populations. Given our assumption that bulges 
form during mergers, these galaxies are bulgeless because they had 
a relatively quiet merger history. Therefore, bulgeless galaxies are 
more likely to form in dark matter haloes hosting few satellites, 
where galaxy merger rates are low enough to assure that the merger 
channel is inefficient in forming a large bulge. The comparison 
between results of the two models, highlights the importance of 
the as sumed prescription for merger times (see also be Lucia et al.l 
[20T0I) . 

We com pare our model predic tions with observational re- 
sults from the lKormendv et ^ ( |2010|) sample. The most interesting 
discrepancy between models and data is found for galaxies with 
10^° < M*/M0 < 10^^, where the fraction of bulgeless galax- 
ies predicted by theoretical models is systematically lower than the 
observational estimate. The discrepancy is significant when com- 
paring predictions from the standard model to the observed frac- 
tion of 'pure discs'. On the other hand, the predictions of the pure 
merger implementations of our SAMs are statistically consistent 
(within 95% confidence level) with the observational estimates, if 
we classify as bulgeless also galaxies with a relevant pseudo-bulge 
(as done in the study by Kormendy et al.). Our results show that 
model galaxies without a classical bulge are not rare: they can ac- 
count for up to 14% of the total mass budget of galaxies in the 



range lO'^^ < M^/M© < 10^^. Our conclusions a re confirmed by 
the qu antitative comparison with recent results by iFisher & Drorvl 
| |20T l|), and by a qualitative compari sons of model predic tions 
with larger galaxy sa mples defined in IWeinzirl et all ( I2OO9I) and 
ICameronet"^ ( l2009l) . These results are in line with previous stud- 
ies in showing the importance of the adopted modelling of disc in- 
stability to account for the estimated B/T of our own Milky Way 
jPe Lucia & Helmill2008l : iMaccio et al1l20ld) . 

We note that our conclusions depend on the assumed B/T 
threshold for the definition bulgeless galaxies. In particular, the 
fraction of bulgeless galaxies decreases when adopting a lower 
threshold (and viceversa). However, for 10^" < M*/M0 and 
M*/M0 > 10^^ our results do not significantly change when 
varying the adopted threshold. This means that, for decreasing 
B/T thresholds, the slope of the fhi{Ah) distribution increases, 
and the nbi(A^sat) distribution becomes more skewed towards low 
Nsat values. Redshift distributions are unaffected. The choice for 
the B/T threshold may thus represent an important aspect of 
our comparison between data and model predictions. Given the 
observational uncertainties in the bulge/disc decomposition (e.g. 
iTasca & Whit3l2005l : iLaurikainen et al II2007I) . we believe that our 
B/T < 0.1 choice is a reasonable one. 

The main result of this paper is that SAMs predicts enough 
merger-quiet galaxies to be (almost) consistent with the dearth of 
classical bulges in the local Universe. In particular, we find that the 
claimed tension between observational estimates based on the local 
volume and predictions from SAMs weakens considerably when (i) 
disc instability is switched off leaving mergers as the only channel 
to create bulges in models, and (ii) only classical bulges are con- 
sidered in the data. Some caveats should be considered. The the- 
oretical models adopted in this paper are tuned to reproduce the 
observed fraction of stellar mass locked in bulges at 2; ~ 0. If 
the instability channel is switched off, these models might under- 
produce the stellar mass in bulges with respect to the observational 
estimates. In addition, as discussed extensively in Paper I, our mod- 
elling of disk instability is very simplistic which has important con- 
sequences for the formation and statistics of bulgefl 

The models considered here do not predict the detailed prop- 
erties of discs (other than baryonic masses and scale-lengths), and 
they assume no angular momentum dissipation, i.e. perfect conser- 
vation of the angular momentum imprinted in the baryonic com- 
ponents (usually modelled assuming some realistic distributions 
derived from numerical simulations). Simulations show that this 
approach is probably oversimplified, and that also haloes with a 
quiet merging history acquire only a fraction of the angular mo- 
mentum required for rotat ional support by tidal torques (see e.g 
iD'Onghia & Burkerdl2004l) . 

Finally, model bulges assemble their mass as the result of both 
mergers and disc instabilities. Given the assumed tight connection 
of classical (pseudo) bulges with mergers (disc instabilities), we 
expect bulges to be composite systems. The KIO sample that we 
have used in this study contains a relevant population (5 out of 
19) of galaxies hosting a relatively large (0.1 < B/T < 0.4) 
pseudo-bulge, and with no evidence of a classical component. In 
order to use these observations to constrain model predictions, 
it is of critical importance to determine to which extent a mas- 



It is worth mentioning that the iBower et alj j200(ih model, which is not 
considered in this work, adopts yet different assumptions for instabilities: 
at each instability event, the whole disc is destroyed and all its baryons are 
given to the spheroidal component. 
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sive pseudo-bu lge may hide a classical component and vice-versa. 
iGadottil lioosl) results suggest that "composite" (actually classical 
bulges with signs of star formation activity) bulges are indeed com- 
mon in the local Universe: unfortunately, the decomposition of a 
galaxy's photometry into the contribution of different components 
is a challenging task, and an unambiguous separation of classical 
and pseudo bulges is curr ently possible only at low redshift (see 
also iTasca & White! I200I) . Advances in this field will enable us 
to increase our knowledge of the complex interplay of the physi- 
cal mechanisms responsible for the distribution of observed galaxy 
morphologies. 
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